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ABSTRACT 

Aims. Drawing an analogy with Active Galactic Nuclei, we investigate the one-zone SSC model of Gamma Ray Bursts 
afterglows in the presence of electron injection and cooling both by synchrotron and SSC losses. 

Methods. We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and 
photons, obtaining the multi-wavelength spectrum as a function of time. We back up our numerical calculations with 
analytical solutions of the equations using various profiles of the magnetic field evolution under certain simplifying 
assumptions. 

Results. We apply the model to the afterglow evolution of GRBs in a uniform density environment and examine the 
impact various parameters have on the multiwavelength spectra. We find that in cases where the electron injection 
and/or the ambient density is high, the losses are dominated by SSC and the solutions depart significantly from the 
ones derived in the synchrotron standard cases. 

Key words, gamma-rays: theory - acceleration of particles - radiation mechanisms: non-thermal 



1. Introduction 

Gamma Ray Bm^st (GRB) afterglows are thought to be 
produced in the Rclativistic Blast Waves (RBW) associ- 
ated with the initial GRB explos ion. Accordin g to the stan- 
dard model (for a review, see iPiranI (120051) 1. the RBW 
expands in the circumburst material and, after sweeping 
some critical amount of mass, it starts decelerating. At 
the same time electrons are postulated to energize at the 
shock front and radia te through synchrotron ( Sari et al.l 
119981: iDermer fc Chia ng 199 ^ or synch r otron and SSC 
(IChian g fc Permed Il999l: iFan et al.l l2008t IWaxmanI [19971: 
Panaite scu fc Kumaill2000f l. thus producing the afterglow 
emission. 

The above picture carries certain analogies to the radi- 
ation models put forward during the last decade to explain 
the multiwavelength (MW) observations of blazar emission, 
and, espe cially, to the 'on e-zon e' SSC model s (for r ecent re- 
views see iBottcheil (|2007t ) and iMastichiadisI (|2009D V These 
models are based on solving an equation for the electron 
distribution function including synchrotron and SSC losses 
and, at the same time, calculate the radiated photon spec- 
trum. This approach allows for time dependency to be taken 
explicitly into account (and thus it can address, for exam- 
ple, blazar flaring), it can treat the non- linear cooling as- 
sociated with SSC and, moreover, it is self-consistent. 

Motivated by these developments, we have applied the 
above technique to the GRB afterglows. Aim of the present 
paper is not to fit spectra or lightcurves, but to focus on 
the impact the various parameters, customarily used by re- 
searchers in the GRB field, have on the MW spectra. We 
present also analytical solutions of the electron equation, 
under certain simplifying assumptions, for a power-law elec- 
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tron injection suffering synchrotron losses. This, amongst 
others, allowed us to test the numerical code. For definitive- 
ness we restrict our analysis to the uniform density case. 

The present paper brings certain improvem ents over 
past efforts in the field (jChiang fc Dermei1ll999r ). First, in 
the analytical solutions we consider both standard and non- 
standard magnetic field evolution. We also use the full emis- 
sivity for synchrotron radiation instead of the (5— function 
approximation. As far as the numerical part is concerned, 
we show that the inclusion of SSC both as electron energy 
loss mechanism and radiation process can bring, under cer- 
tain circumstances, significant departures from the stan- 
dard solutions which include only synchrotron radiation as 
an energy loss mechanism. Finally, we have included for the 
first time, to the best of our knowledge, 77 absorption both 
as a 7— ray attenuation and particle reinjection mechanism. 

The paper is structured as follows. In §2 we review the 
basic hydrodynamics and radiative concepts of the GRB 
afterglow adopted by the extensive literature on the sub- 
ject. We also make an analogy to the blazar case. In §3 
we present the analytical results for various profiles of the 
magnetic field. In §4 we present the numerical code and 
the tests we have performed to check its validity. In §5 we 
show various numerical results and we conclude in §6 with 
a summary and a discussion. 

2. Physics of the Relativistic Blast Wave 

2.1. Hydrodynamics 

We assume a shell of material with initial mass AIq moving 
with initial bulk Lorentz factor Fg. This will sweep am- 
bient matter and will start decelerating with a rate de- 
termined by energy and momentum conservation. The re- 
sulting blast wave is modeled as having a cross sectional 
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area, A{r), that depends on the distance (measured in the 
frame of the explosion) . In our work we consider a spherical 
blast wave with A{r) — Anr'^ sweeping a constant density 
ambient matter. We also assume that the bulk kinetic en- 
ergy which is converted to internal energy is not radiated 
away, contributing to the inertia of the blast wave (non- 
radiative limit). In this case the deceleration of the blast 
wave is determined by a pair of ordinary differential equa- 
tions (jBlandford fc McKedfTOTj) : 



dT 
dr 
dM 
dr 



A(r)p(r)(r^(r)-1) 
M(r) 

pir)A{r)T{r) , 



(1) 
(2) 



where F — r(r) is the bulk Lorentz factor of the mate- 
rial, M = M{r) is the total mass including internal kinetic 
energy and p — p{r) is the mass density of the ambient 
matter. Note that Eqn. (2) implies a non-radiative blast- 
wave evolution. We will adopt this assumption throughout 
the paper. 

It can be shown (jBlandford fc McKedll976t) that the 
above system has an analytic solution given by: 



r(r) = 



A(.t3 - 1) + To 



(1 + 2roA(a;3 - 1) -f- ^{x^ - l)2)i/2 



where A 



^ and Mo = 



(3) 



Here rg is the 



3Mo ' - ro 

initial radius of the blast wave and po is the mass density 
of the uniform ambient matter. The evolution of the bulk 
Lorentz factor can be separated into three regimes. The first 
regime corresponds to the initial period of free expansion of 
the blast wave, during which T[r) « Fq. The decelerating 
phase follows, where the bulk Lorentz factor can be modeled 
as 



r(r) 



where 



ro 

2 



Rd 



-3/2 



10^3/2 

2 



and 



Rd 



3Eo 



\47rnmpc2rQ 



1/3 



(4) 



(5) 



(6) 



is the decelera t ion radius of the blast wave 
(jRees fc Meszarod Il992| l. This power-law r'^/^ 
pendence of the bulk Lorentz factor is often quoted for 
a non-radiative blast wave decelerating in a uniform 
medium. During the third and final regime the blast wave 
is practically non-relativistic. 

Energy conservation gives the rate of accreted kinetic 
energy in the lab frame: 



dE 

— = c^A{r)p{r)P{T^{r)^T{r)) 



(7) 



This expression follows from the equation of motion (eqn. 
([1])) and is applied regardless of whether the blast wave is 
in the radiative or adiabatic regime. As the rate of energy 
accreted is a Lorentz invariant, expression ^ holds also in 
the comoving frame of the blast wave. 



2.2. Radiation 

The RBW does not only sweep matter and decelerates, as 
was discussed above, but it is assumed to be able to ener- 
gize particles as well. While there are no detailed models as 
yet to explain the way particles achieve high energies, it is 
assumed that an ad-hoc fraction Ce of the accreted kinetic 
energy is injected into non-thermal electrons with a power- 
law form. A second assumption concerns the lower and up- 
per cutoffs of the electron distribution 7min and 7niax- Since 
the normalization go of the electron injection is set by the 
relation 



907 ^(7 - l)meC^d-f 



dE 
'~dt' 



(8) 



where 7 is the Lorentz factor of the electrons and p is the 
electron injection spectral index, it is evident that for p > 2 
(as is typically assumed) , only the lower limit 7min will play 
a role in determining the integral in Eq.® and thus go- 
Therefore the choice of 7inin has important consequences 
on the results. This can be at least of the order of F in 
the comoving frame (jMastichiadis fc Kazanasll20 09') but in 
case where the electrons are in equipartition with protons 



it can be a factor of ( mp/we) higher, i.e. 7n 



(rnp/me)r 



(jKatz fc Piranlll997t|Panaitescu fc Meszaroslll998[ ). 

Since the energetic electrons will emit synchrotron ra- 
diation, a prescription for the magnetic field is also re- 
quired. This again is a source of major uncertainty. The 
usual assumption is that the magnetic field is in some type 
of equipartition with the particles; this implies that the 
magnetic energy density takes a fraction cb of the mass ac- 
cumulated in the RBW, so the magnetic field is given by 
the relation 



B = (327rnTOpeBc2)^/^F. 



(9) 



However one can consider differen t types of behavior, like 
B ~ r-i (|Vlahakis fc Konig]|l2003l) . 

Depending on the value of the magnetic field, the elec- 
trons can cool (i.e. radiate all their energy) in a dynami- 
cal timescale or remain uncooled, (i.e. keep their energy). 
These two cases have been called 'fast' and 'slow' cooling 
(jSari et al.lll998D . From the standard sol utions of the elec- 
tron kinetic equations (|Kardashevlll962l ) it is known that 
cooled electrons have steeper distribution functions than 
uncooled ones, i.e. cooled electrons have an energy depen- 
dance that is proportional to while uncooled ones 

have a 7"^ dependance, i.e. they still retain the spectrum 
at injection. 

High energy electrons can also lose energy by inverse 
Compton scattering on ambient photons. These photons 
can illuminate the source of electrons externally or they can 
be the synchrotron photons mentioned earlier. In this lat- 
ter case the process is called synchro-self Compton (SSC). 
The inclusion of the inverse Compton process has two im- 
portant results: (a) It will produce a high energy spectral 
component in the photon spectrum and (b) depending on 
the respective magnetic and photon field energy densities, 
it can alter the electron distribution function, thus affecting 
directly the shape of the radiated photon spectrum. 



2.3. Comparisons to blazar models 

The radiation coming from the GRB shock, as was de- 
scribed above, has the same underlying physical principles 
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with the so-called 'one-zone' SSC AGN leptonic models 
set f orward t o explain the MW spectrum of blazar emis- 
sion (jlnoue fc Takah ara 1996; Mastichi adis & Kirk ,1997) 
These mod els, based on earlier i deas s et bv iMaraschi et afl 
(I1992D andU oom fc Marscherj (|1996( ). address essentially 
the same problem of electron injection, cooling and photon 
radiation. However they do not take a 'ready' electron dis- 
tribution but, rather, they obtain it from the solution of a 
kinetic equation which contains injection, radiative losses, 
physical escape from the source and possibly reinjection 
of particles as secondaries from photon-photon absorption. 
The electron equation is coupled to an equation for pho- 
tons that has the usual synchrotron and inverse Compton 
emissivities written in a way as to match the radiative elec- 
tron losses. It can also have extra terms such as synchrotron 
self-absorption, photon-photon pair production, etc. 

This kinetic equation approach has the advantage that 
it is self-consistent, i.e. the power lost by the electrons is 
radiated by the photons. Moreover the photon energy den- 
sity can be calculated at each instant and this feedbacks 
though the SSC losses on the electron equation, thus this 
approach can treat the SSC intrinsic non-linearity. 

Therefore it would have been instructive for one to con- 
struct a similar approach to model the radiation of GRB af- 
terglows. However, one should have in mind, that in GRBs, 
despite the physical analogies to AGNs, the situation has 
some obvious differences which come mainly from the hy- 
drodynamics of the GRB outflows as this was outlined in 
2.1. 

1. AGN modeling involves usually stationary states. 
This means that the the radius of the source Rs , the Lorentz 
factor r and the magnetic field strength B are all consid- 
ered constant. This holds even when short flares are mod- 
eled in a time-dependen t wav ((Ma stichiad is fc Kirld 119971 : 
iKrawczvnski et al.l [200l iKatarzvnski et al.l \20oW as we 
mentioned earlier, in GRBs all of the above are functions 
of the distance r from the origin of the explosion. 

2. In GRBs once the profile of the external density is set, 
then the injection of electrons has at least an upper limit 
as Ce cannot exceed unity. There is no such constraint for 
AGNs as the injected power is essentially a free parameter. 

These differences imply that the numerical codes devel- 
oped for blazars cannot be used as they are for the GRB 
afterglows, but they have to be modified to take the above 
into account. We will present such a code in §4. However we 
first show some analytical solutions that we have derived 
for the coupled hydrodynamic-radiation problem. 



3. Analytical solutions 

In this section we will first present the kinetic equation for 
the electron distribution function, then we will show the so- 
lutions for different magnetic field configurations assuming 
that synchrotron losses dominate and we will end by using 
these solutions to derive the slope of the synchrotron and 
SSC lightcurves which corresponds to each B-field configu- 
ration considered. 



3.1. Kinetic equation of electrons 

The equation which governs the electron distribution is: 



(10) 



where iV(7, t) is the number of electrons having Lorenz fac- 
tors between 7, 7 -I- fi7 at a time t as measured in the co- 
moving frame. The second term of the left hand side takes 
into account energy losses of the relativistic electrons due 
to synchrotron emission and inverse compton scattering in 
general. The term on the right hand side of the equation 
describes electron injection. Equation (fTO|) can be also ex- 
pressed in terms of the distance r measured in the frame of 
the explosion, instead of t, through dt 



dr 



dNij,r) 
dr 



-l|-(5(7,r)Ar(^,r)) = | 



Qh,r) (11) 



where iV(7, r) is now the number of electrons having Lorenz 
factors between 7, 7 -f ^7 at a radius r as measured in the 
frame of the explosion. In order to obtain some analytical 
results of the above equation we consider only synchrotron 
losses. We assume that a fraction ec of the accreted kinetic 
energy is injected into non thermal electrons. The injection 
of the relativistic electrons can be modeled as a power-law 
with arbitrary choices of minimum and maximum Lorentz 
factors in energy: 

Q(7, r) = go(r)7"^e(7 - 7min)e(7„,ax - iMr - tq). (12) 

At each radius the normalization of the electron distribu- 
tion is given by the prescribed fraction of the power avail- 
able as bulk kinetic energy. Equation ([5]) leads to 



90 (0 



(r^ - ry 



(13) 



The constant which appears in equation (|13p is given by: 

ifp^2 



47rpoc 



^ ( "'min 
1 p-2 



^ mcC^ ( 1: 



+ 



lip ^2 



3.1.1. Magnetic Field of the form B{r) = Bor-^/'^ 

As was briefly discussed in Section 2, this type of B-field is 
the one customarily adopted for the GRB afterglows. Eqn. 
([9|) implies that a fraction ee of the accumulated mass on 
the RBW goes to amplify the B-field. As F c>c r^-^/^ in the 
decelerating phase (c.f. eq.(I3])) the above prescription for 
the B-field is derived at once. Thus Eqn.® can be further 
written 



B = Bor 



'3/2 



with 

Bo ^ (327rnmpeBC^)^^^f 



(14) 



(15) 



where the constant F is defined in equation Then the 
term of synchrotron losses becomes 



6(7, r) 



-ao- 



7 



where an — „ 
simple expression in case oi p — 2 

2k f 



(16) 

The solution of equation pT|) has a 



(17) 
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expressions for the electron distribution in these regimes. 




logy 



Fig. 1. Comparison between analytical (solid line) and nu- 
merical (dotted line) solution in case of p = 2 at radius 
r = 1.5 X lO^^cm or at observer time 2.7 x 10^ s. There 
is practically no difference between the two solutions. The 
parameters used for this plot are: = 0.1, eb = 0.001, n = 



Icm-^To = 400, £;o = 10^-^erg s-\-f, 
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where 

K{'y,r) 

and 

7b = 



„3/2 _ 3/2 



if 7 < 7b 



„3/2 _ 



1 


f 2ao 




7max 


cT 





(18) 



We point out that the second branch of the above solution 
does not necessarily describe a cooled electron distribution. 
Only if the relation 



2ao7 



< 1 



(19) 



holds, then the electron distribution can be considered 
cooled, i.e. N oc 7"^. A more physical approach is 
to consider an electron with Lorentz factor 7c, which 
cools in a timescale equal to the dynamical ^dyn 
(jPanaitescu fc: Kumaill2000l ). This is given by: 



2j_ 

5 cT 



7c 



(20) 



Then, electrons with Lorentz factors greater than 7c cool 
sufficiently. Condition which arose mathematically, 

is equivalent to 7 7c. This condition and relation (fT9|) 
differ only by a factor of 5. From the above discussion it is 
clear that the Lorentz factor 7b must not be confused with 
7c, that determines at which Lorentz factor the cooling 
of the distribution becomes dominant. An example is 
shown in Figure (TJ where the change of the slope in the 
power law begins at 7 « 5 x lO'*. For the parameters used 
7c = 5.2 X 10* and 7b = 2.8 x 10^. In the case where p ^ 2 
the integral Ip can only be estimated in two asymptotic 
regimes. The complete calculations and approximations 
can be found in the Appendix. Here we present the 



In the uncooled regime the relation 7 <JC 7c holds 
and the solution takes the form: 



Died 



2A:cr 
3c 



(21) 



As we mentioned above in the cooled regime the condition 
7 ^ 7c holds. Then 



krT^ 



□led 



aoip- 1) 



(22) 



3.1.2. Magnetic field of the form B{r) = Bq 



This type of magnetic field might also be rel ated to 
GRB outflows (see e.g. IVlahakis fc Konigll (|2003f )). Here 
we present the analytical solutions of equation PT|) for this 
type of magnetic field, while the complete calculations can 
be found in the Appendix. Again if p = 2 the solution has 
the simple form 



N{j,r) 
where 

and 

7b = —i 



3c7^ 



(23) 




if 7 < 7b 



7b 



1 



The constant ao is given by 
arBirl 



(24) 



(25) 



If p > 2 then the solution is found in two regimes , as in 
(13X11) . 

The uncooled part of the electron distribution is given by: 
2fccf 



uncooled ' 



3c 



(26) 



which is exactly the same as the one calculated for the 
magnetic field B cx r^^/^. In this regime the cooling 
timescale of electrons is much greater than the dynamical 
timescale. Thus, this part of the electron distribution will 
not be affected by a different type of magnetic field which 
is the cause of electron cooling. 

The situation is different for the part of the distribu- 
tion where cooling is dominant: 



N, 



cooled 



ao{p ~ 1) 



r7 



(27) 
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3.1.3. Constant magnetic field Bq 

In order to check our analytical results we have also solved 
equation (jlip for the case of constant magnetic field. This 
calcu lation has already been done by iDe rmer fc Chiang 
(|1998|) and thus we can compare our results with theirs. 
The outline of the comparison can be found in the Appendix 
l|A.3|) . In the uncooled regime the solution is given by: 



Table 1. Negative of slopes of Time Profiles 



uncooled 



3c 



(28) 



The above expression for the uncooled part of the electron 
distribution is again the same as for the other types of mag- 
netic field presented in the previous sections. In the cooled 
regime 



N, 



cooled 



aoip- 1) 



(29) 



In this case the cooled part of the distribution reduces as 
the radius of the blast wave increases. This behavior of the 
cooled part of the distribution differs from the one pre- 
sented in the previous sections where the total number of 
electrons within the shell increased with increasing radius. 

3.2. Analytic flux time profiles 

The kinetic equation of the electron distribution is being 
solved in the comoving frame, as shown in section (jS.ip . 
Synchrotron and SSC spectra are also first calculated in 
the comoving frame and then transformed into the observer 
frame. For this we use a relation which connects time in 
the observer frame and radius r which appears in all our 
analytical solutions. Thus, 



^obs 



dr 
2d^ 



(30) 



where r is the radius of the blast wave measured in the 
comoving frame. If the distance of the source from the ob- 
server is D then the respective synchrotron and SSC fluxes 
at the observer frame are given by: 



-fs(^obs! ^obs) 



<-^s(^obsi ^obs)r 



(31) 





B — const 


B OC r"^ 


Boc r-'^'l'' 


uncooled 


Sp-3 
16 


5p-l 
16 


3p 
8 


cooled 


3p+10 
16 


5p+4 
16 


6p+l 
16 


uncooled 
^ssc 


3p-l 
16 


5p+l 
16 


3p+l 
8 


cooled 
L^ssc 


3p+22 
16 


5p+8 
16 


6p+l 
16 



where Ce — Cc{r) is the normalization factor of the elec- 
tron distribution. After the deceleration radius (eqE]) it 
is straightforward to show that tobs oc r"*. Thus the syn- 
chrotron and SSC fluxes can be expressed as power laws 
of the observed time with exponents agjCtssc respectively. 
The exponents are found in the asymptotic regimes of un- 
cooled and cooled electron distribution and are presented 
in Table [TJ The slopes of the synchrotron flux time pro- 
files in the non-radiative regime (see Table [J) coincide wit h 
the respective ones presented bv lDermer fc Chiangj (jl998[ ). 
Table [T] shows that SSC flux time profiles are steeper than 
the respective synchrotron profiles both in the uncooled 
and cooled regime, for all the magnetic field configurations 
discussed in section [XT] 



4. Numerical approach 

4.1. The code 

The equation solved in the previous section forms the ba- 
sis of the electron kinetic equation; we proceed now to 
augment this with more processes and to solve it numeri- 
cally. Since SSC losses depend on the synchrotron photons 
energy density, we have to write an accompanying equa- 
tion for photons which is coupled to the electron equa- 
tion. S imilar type of equations have been solved for the 
blazar llMastichiadis fc Kirk 1997) and prompt / early after- 
glow GRB cases (jMastichiadis fc Kazanasl2009l) - note that 
in this latter case a third equation for protons was added. 

Assuming, as before, that the electrons are a function of 
distance from the center of the explosion and energy, their 
equation reads 



p / J. \ _ '^ssc(^^obsi ^obs)r , , 

where Jg, Jgsc are the synchrotron and SSC power per unit 
frequency emitted in the comoving frame. It is interesting 
to examine the dependence of the observed fiuxes on time, 
in cases where the electron kinetic equation can be analyt- 
ically solved. Thus, we consider the cases which we have 
treated in (|3.ip . Since we are interested only in the time 
dependency, in what follows we can work using proportion- 
alities. 

Equations and (|32p can be reduced to the simplified 
form 

F, OC C,T^+°'B'^ (33) 
Fssc OC ^ri+"S^ (34) 



+ CT' + £r + = Q^e' + QV (35) 
while the corresponding photon equation is 



^ + ^ + + £7 = Q^y" + Qi;^ (36) 

OT tcr 

The operators C denote losses and escape from the sys- 
tem while Q denote injection and source terms. The un- 
known functions rie are n-y are the differential number den- 
sities of electrons and photons respectively and the phys- 
ical processes which are included in the kinetic equations 
are: (1) electron synchrotron radiation and synchrotron self 
absorption (denoted by the superscripts "syn" and "ssa" 
respectively); (2) inverse Compton scattering ("ics"); (3) 
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photon-photon pair production ("77") and (4) adiabatic 
losses ("adi"). 

The numerical code keeps the same philosophy, as far 
as the phy sical processes are concerne d, with the one de- 
scribed in iMastichiadis fc KirH (|1995D (hereafter MK95). 
However since various modifications have been introduced, 
we summarize briefly the expressions used: 

(i) Synchrotron Radiation: (a) The electron loss term 
£1^" is given by expression (34) of MK95. (b) The pho- 
ton emissivity term Q^^" is given using the full high en- 
ergy emissivity term (see, e.g. iBlumenthal &: Gouldl [19701 ) 
instead of the delta- function approximation used in MK95. 

(a) Synchrotron-Self Absorption: (a) The photon ab- 
sorption term is used as in MK95 (Eqn. 39). (b) There 
is no matching term for electron heating due to this pro- 
cess. However, as synchrotron self absorption is expected to 
be minimal for the parameters which are of interest here, 
the error introduced by this omission is expected to be neg- 
ligible. 

(Hi) Inverse Compton scattering: (a) The electron 
loss term C^^ is given by solving Eqn. (5.7) of 
IBlumenthal fc Gouldl ([l970l) . (b) The photon emissivity 
term Q^^^ uses relation (2.48) of the same paper. 

(iv) Photon-photon pair production: (a) The electron 
injection term QJ'' is given by expression (57) of MK95. 
(b) The photon absorption term is given by expression 
(54) of MK95. 

(v) Electron injection : The quantity Q"^' is the electron 
injection rate which can take any functional form of dis- 
tance r and energy 7. Following the usually assumed case, 
we take it to be of a power-law form as given by Eqn ([T^ . 
The power-law index, the normalization and the upper and 
lower energy cutoffs can be treated as free parameters. 

(vi) Photon escape: This is characterized by the light 
crossing time of the source tcr = Rs/c in the comoving 
frame. This is not constant but changes according to the 
relation Rs = r/T. 

There is still one more free parameter to be addressed 
for specifying the code quantities and this is the behavior of 
the magnetic field with radius. While it is trivial to adopt 
any r-dependence we will use the standard case B cx r~^/^ 
except in one of the tests that follow. 



4.2. Tests 

There are various tests that we have performed to test 
the validity of the code. As the various rates of the ra- 
diative processes have essentially the same form as the 
ones used in the past to model AGN MW em ission 
(jMastichiadis fc KirkI 119971: iKonopelko et all l2003l ). they 
have been checked many times a gains t the results of e.g. 
ICoppil (|1992[ ): iKatarzvnski et al.l (|2005[ ) and others. 

The new aspect introduced in the code is its dynamical 
behavior and in order to test this we have checked the code 
extensively against the analytical results (spectral shape 
and lightcurve slopes) derived in the previous section. As an 
example we show in Figure [5] such a comparison of the mul- 
tiwavelength spectra derived with the analytical method of 
the previous section with the results of the numerical code. 

The other major test was to compare with results al- 
ready published in the literature, such as the ones given in 
(|Fan et al.ll2008h . Figure [3] reproduces two curves of Fig. 4b 
of the aforementioned paper with very good agreement. 




log V (Hz) 



Fig. 2. Synchrotron and SSC spectra for a case with Fg = 
100, Eo = IG^^crg s-\ Rq = IQi^cm, n = 1 cm'^, p = 2, 
7min = 10, 7max = 10* and B{r) = 106(fi) Gauss. The 
spectra are calculated at ?■ = 2 x 10^^ cm. The source was 
assumed at a distance of _D = 3 Gpc. Analytical results 
are depicted with full line, while the numerical ones with 
the dotted one. The difference in the low part of the spec- 
trum arises from the fact that in the analytical results we 
did not solve for the electron distribution below 7inin- In 
the numerical case we did not impose such restriction. For 
comparison reasons we ignored SSC losses in the numerical 
code. 




log V (Hz) 



Fig. 3. Multiwavelength spectra with parameters similar to 
the ones stated in Fig. 4b of Fan et al. (2008). These should 
be compared with the two top curves of the aforementioned 
Figure. 



5. Results 

In what follows we apply the numerical code to examine the 
dynamical evolution of the MW GRB afterglow spectra. For 
this we use the standard afterglow adopted parameters, i.e. 
the magnetic field is given by rel. ^ with ee a free param- 
eter, while the injected electron luminosity is controlled by 
Eo which is also treated as a free parameter. Furthermore 
we adopt a value for 7min that is not constant but is of the 
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log V (Hz) 

Fig. 4. Photon spectrum obtained at radius R — 3.4 x 
10^^ cm or equivalently at observer time tobs = 7 x 10"^ s 
for To = 400, Eq = IQ^^eig s-\ n = 1 cm-^ , = .1, 
ee = .001. The electrons were assumed to have a power- 
law distribution with slope p = 2.3 while their maximum 
cutoff was 7max = 4 X lO*". 7inin was given by expression 
(|57|) . The full line corresponds to the spectrum when all 
processes are included, the short-dashed one when 77 is 
omitted and the dotted line one when only synchrotron is 
included. 




Fig. 5. Electron distribution function for p — 2.3, normal- 
ized by the factor aTrricC^ Rg and multiplied by 7^, that 
corresponds to the photon spectrum of Figure HI The full 
line corresponds to the solution when all processes are in- 
cluded, while the dotted line one when only synchrotron 
losses are included. 

form (|Sari et al.lU"998[) 

p-2 nip 

7min = -Co r (37) 

p — 1 rUe 

Our aim is to see whether (i) SSC losses can modify the 
electron distribution function, and therefore, the photon 
spectrum and (ii) photon-photon absorption, a process that 
has been neglected thus far, can be of some importance, not 
only taken as a 7— ray absorption mechanism but also as 
an electron (and positron) reinjection one. 

As a first case we show an example for typical values 
assumed usually for GRB afterglows. Figure [3] shows the 



photon spectrum obtained at radius R — 3.4 x 10^^ cm for 
Fo = 400, n = 1 cm~3, Eo = .1, cb = .001. The electrons 
were assumed to have a power-law distribution with slope 
p = 2.3 while their maximum cutoff was taken to be con- 
stant and equal to 7max = 4 x 10''. Here and in the next 
Figures the GRB was set at z—1. Note also that we assume 
that the evolution of the Lorentz factor F follows the adi- 
abatic prescription implied by Eqn. (2) and that we have 
not taken into account any attenuation for TeV 7— rays due 
to absorption on the IR background. The full line curve de- 
picts the photon spectrum when all processes are included, 
the dashed line one when photon-photon absorption is left 
out and the dotted line one when inverse Compton scatter- 
ing is also omitted both as emission in the photon equation 
and as loss mechanism in the electron equation, i.e. this case 
can be considered as pure synchrotron. Although this latter 
case is clearly an oversimplification, we have included it for 
comparison. One can see that 77 absorption influences only 
the highest part of the spectrum by making it steeper. Pair 
reinjection does not alter significantly the lower spectral 
parts because only a very small fraction of the energy has 
been absorbed and is thus available for redistribution. On 
the other hand, SSC losses have an impact on the spectrum. 
This effect can be seen better in Figure [5] which shows the 
electron distribution function with (full lines) and without 
(dashed lines) SSC losses at the aforementioned radius. As 
in standard theory, the braking energy 7c divides the cooled 
(7 > 7c) from the uncooled (7 < 7c) part of the electrons. 
A first comment one could make is that the synchrotron 
break does not appear as a sharp turnover but as a gradual 
one which affects the power law index of electrons at least 
one order of magnitude around 7c. A second comment is 
that inclusion of SSC losses changes the 'pure' synchrotron 
picture. The electron distribution function becomes flatter 
and this is a result of the SSC losses. The specific shape 
can be explained because the SSC losses at each electron 
energy consist of both losses in the Thomson and the Klein 
Nishina regime. As the electron energy increases, the frac- 
tion of Klein Nishina to Thomson losses also increases, with 
the result the total SSC losses to be reduced. 

5.1. Dynamical Evolution 

Figure ini shows snapshots of photon spectra obtained for the 
same parameters as above at three different radii; R — i?dcc, 
R — 3.2i?dec and R = 10-Rdec- As before, full Hue curves 
depict the spectrum when all processes are included, dashed 
line ones when photon-photon absorption is left out and 
dotted line ones when inverse Compton scattering is also 
omitted. Obviously 77 absorption affects only the highest 
energies and does not play any significant role throughout 
the evolution. 

On the other hand, SSC losses seem to play a role that 
becomes slightly more important as the radius increases. 
This can be understood from the fact that as the magnetic 
field drops, two contradicting results occur: One is that 
synchrotron cooling becomes less efficient (it moves from 
'fast' to 'slow') and the number of available soft photons 
for upscattering is reduced. However, these same photons 
become softer with radius (since both B and 7min are re- 
duced outwards) and therefore one expects more collisions 
in the Thomson regime where electron losses become more 
efficient. This effect can be seen by comparing the shapes 
of the MW spectra at the three radii: the shape of the 
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Fig. 6. Multiwavelength GRB afterglow spectra for Fq = 
400, = IG^^erg s-\ n = 1 cm,-^ , = .1, eb = -001, at 
R = Rdec, R = 3.2i?doc and R = 10i?dcc (top to bottom) 
or at equivalent observer times tobs — 7.2 s, 2.7 x 10^ s, 
2.4 X lO'* s respectively. The electrons were assumed to 
have a power-law distribution with slope p — 2.3 and 
7max = 4 X 10''. Full line curves depict the spectrum when 
all processes are included, dashed line ones when photon- 
photon absorption is left out while dotted line ones are pure 
synchrotron cases. 
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Fig. 7. Multiwavelength GRB afterglow spectra for Fq — 
400, Eo = 1053gj.g g-i^ ^ ^ I ^^-3^ ^ 001. The elec- 
trons were assumed to have a power-law distribution with 
slope p = 2.3 and 7max = 4 x 10*". All photon spectra are 
calculated at the radius R = 3.2i?(joc or at observer time 
t^y^^ = 2.7 X 10^ s with ee = 0.1, 0.01 and 0.001 (top to bot- 
tom) . Full line curves depict the spectra when all processes 
are included and dotted line ones depict the correspond- 
ing spectra when only synchrotron radiation is taken into 
account. 

-9 I 1 1 1 1 1 1 1 1 1 1 



SSC component starts resembling the synchrotron one as 
the distance increases and the SSC losses are dominated by 
collisions in the Thomson regime. 

5.2. Role of te 

The effect that has on the spectrum is more straight- 
forward. As Ec decreases (for fixed ee) the electron spec- 
tra are increasingly dominated by synchrotron losses and 
the effect of SSC losses becomes marginal. This is shown 
in Figure [7] which depicts the radiated photon spectra for 
three values of = 0.1, 0.01 and 0.001 (top to bottom). 
Here ee = 0.001, while the rest of the parameters are as 
in the previous case. In order to avoid confusion we have 
calculated all spectra at radius R — 3.2i?dcc- Note that, as 
Ce decreases, the SSC component drops as the quadratic of 
the synchrotron component, a fact that is well known in 
the SSC AGN models. 



5.3. Role ofcB 

Figure [5] shows the effects that eb has on the MW spectra. 
Here the run has the same parameters as before, however 
the spectra are calculated at the same radius R = 3.2_Rdoc 
with eb =0.1, 10~^ and 10~^. SSC losses change the spec- 
tra only for intermediate values of cb- High values of eb lead 
to fast cooling and a predominance of synchrotron radia- 
tion. On the other hand, low values of eB lead to inefficient 
cooling and only the high synchrotron frequencies are af- 
fected by the SSC cooling which occurs in the deep KN 
regime. 



-10 - 
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Fig. 8. Multiwavelength GRB afterglow spectra for Fq = 
400, Eo = 10^3gj.g g-i, n = 1 cTO-3, = .1. The elec- 
trons were assumed to have a power-law distribution with 
slope p = 2.3 and 7inax = 4 x 10^. All photon spectra are 
calculated at the radius R = 3.2i?dcc or at observer time 
tohs = 2.7 X 10^ s with EB = 0.1, 10"^ and lO"'^ (^p to 
bottom). Full line curves depict the spectra when all pro- 
cesses are included and dotted line the corresponding spec- 
tra when only synchrotron radiation is taken into account. 

5.4. Role of external density n 

Figure [5] shows the way the MW spectra change in the case 
when the density is increased to 1000 part/cm^. This figure 
has to be directly compared to Fig. ID Inclusion of the SSC 
losses makes the spectrum to depart significantly from the 
pure synchrotron case. Therefore this is a clear Compton 
dominated case with the SSC component exceeding the syn- 
chrotron one by an order of magnitude. Furthermore 77 
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Fig. 9. Multiwavelength GRB afterglow spectra for Fq = 
400, Eo = 1053gj.g g-i^ „ ^ iQQQ ^^-3^ ^ I j^j^j 

cb = 0.001. The electrons were assumed to have a power- 
law distribution with slope p = 2.3 and 7max = 4x10^. All 
photon spectra are calculated at the radius R = 2.7 x 10^^ 
cm or equivalently at observer time iobs = 2.8 x 10^ s. The 
full line curve depict the spectrum when all processes are 
included, the dashed line one when photon-photon absorp- 
tion is left out while the dotted line one is a pure syn- 
chrotron case. Notice that the inclusion of the SSC losses 
redistributes the radiated power to high energies while 77 
has an effect on the entire spectrum through pair injection. 

absorption produces a contribution that affects the entire 
spectrum as substantial pair injection redistributes the lu- 
minosity from the high energy end to lower. 

6. Summary/Discussion 

In the present paper we have applied the 'one-zone' SSC 
model which is customarily applied to the multiwavelength 
blazar emission to the GRB afterglows. For this we have 
used a numerical code that treats electron injection and 
cooling and we have calculated self-consistently the elec- 
tron distribution and radiated photon spectrum at each 
radius of the relativistic blast wave. In this sense the 
present work should be cons idered as complimentary to 
IChiang fc Permed (|1999f ) and lFan et all (|2008). 

One difference between the modeling of blazar emission 
(even in flaring conditions) and GRB afterglows is that in 
the latter case there is continuous evolution of the bulk rela- 
tivistic Lorentz factor F, the magnetic field strength B and 
the radius of the emitting source i?s = r/F. In order to test 
the numerical code in this new setting we have solved an- 
alytically the kinetic equation of the evolving electrons for 
the uniform density case under certain simplifying assump- 
tions and compared the resulting electron distribution func- 
tion and radiated spectra to the ones given from the code 
for similar parameters. The analytical solutions, details of 
which can be found in the Appendix, have their own inter- 
est despite the introduced simplifications as they can give 
simple expressions for the dependence of the lightcurve on 
the spectral index in the case where the losses come solely 
from synchrotron. We find, for example, that the magnetic 
field prescription plays a significant role in determining the 
slopes of the flux time profiles (see Table 1). Moreover they 



treat correctly the synchrotron cooling break which should 
not be taken as an abrupt change of slope but as a very 
gradual one. 

Restricting ourselves to the uniform density case in the 
adiabatic approximation, we have also performed numerical 
runs of the standard GRB afterglow approach as this can 
be summarized by the introduction of the usual parameters 
cb and e^- Our aim was two- fold: (i) To see whether GRB 
afterglows can run into Compton dominant cases similar 
to the ones found in modeling AGNs and (ii) to investi- 
gate the effect 77 absorption and subsequent pair reinjec- 
tion might have on the afterglow multiwavelength spectra. 
Starting from the latter case first we found that, at least for 
the set of parameters used, 77 absorption has an effect, as 
a photon attenuation mechanism, only at the very high en- 
ergy part of the spectrum. Therefore it could affect only the 
potential TeV GRB observations - we point out also that 
TeV emission from GRBs at high z will be subject not only 
to internal 77 attenuation, but also on extra absorption by 
photons of the IR background, an effect that already has 
important consequences even for low z TeV blazars (see, 
e.g. iKonopelko et all (|2003t )) However, as a pair injection 
mechanism 77 absorption has an effect over the whole spec- 
trum that becomes increasingly important with increasing 
values of the external density (compare Fig[3]with FiglH]). 

At the same time, we found that SSC has an impact on 
the MW spectra as it introduces an extra source of cooling 
on the electron distribution. As one would expect it depends 
also on the parameters Cc, eb and, more critically, on the 
value of ambient density n. For low values of Eo it does not 
play practically any role and the cooling comes solely from 
synchrotron radiation (see Fig[7|). Thus the GRB afterglow 
is in the "synchrotron dominated" regime. However for high 
values of Co it becomes important and this can be seen from 
the fact that in the vF^, spectra the SSC component carries 
about the same luminosity as the synchrotron component 
even in the n — 1 cm~^ case. The effect becomes quite 
severe for higher values of n where the SSC component 
carries an increasingly higher part of the luminosity (see 
Fig. [5]). These afterglows are clearly Compton dominated. 

The fact that we have limited our analysis to the adia- 
batic case restricts the allowed values of Cg- As a limiting 
case we have adopted the value Co = 0.1 which introduces 
a maximum (corresponding to the fast cooling case) error 
of a few percent in the Loren tz factor F as this has t o be 
corrected for radiative losses (jChiang fc Dermei1ll999[ ). In 
any case, these corrections are small, they do not affect the 
spectral shapes and, therefore, they cannot alter the basic 
results of the present paper. 

Concluding we can say that the inclusion of the effects of 
SSC (not only as an emission, but as a loss process as well) 
and, to a smaller degree, of 77 absorption can have a signif- 
icant impact on the GRB afterglow multiwavelength spec- 
tra. As expected, the overall picture that emanates from 
these processes is rather complicated. However, as a rule of 
thumb, we can say that the aforementioned processes be- 
come increasingly important for high external densities and 
high ec cases. Their inclusion has some non-trivial conse- 
quences in the lightcurves as both the photon spectra slope 
and the energy break are modified. On the other hand, they 
can be rather safely neglected for high eb and low Ec or n 
cases. 

Based on the above results we expect that SSC losses 
will play an important role in a wind-type density profile. 
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especially at early times when densities are high. This will as mentioned in section ()3.1.1|) . The integral Ip is then 
be the subject of another paper. 



calculated 
1 
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(A.IO) 



Appendix A: Analytical solutions of the electron 
kinetic equation 

A.l. Magnetic field of the form B = Bor"^/^ 

The solution of equation (llip is 
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or changing the integration variable 
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7V(7,r) 



Q:o7 



The electron distribution is in general given by: 



7V(7,r) 
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(A.12) 




where constant ao is defined in section [3.1.1l 



(A.13) 



where F(a, b; c; z) is the Gaussian hypergeometric function. 
The hypergeometric function has a power series represen- 
tation. The three first terms of the series are: 
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The condition << 1 allows us to expand each term 

We define ^ = The integral for p > 2 (as it is typically of equation(|X9| in terms of We work only in zeroth 

order, in order to obtain the simplest expression for Ip and 
its depcndance on r, 7: 



assumed) can be estimated in two regimes: 
Uncooled 

In this regime 7 < 7c and 7, < 7max- Taking also 
into account that 7* > 7 and p > 2 equation (|A.4p 
becomes: 
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The term of equation (|A.17p will not be taken into account 
later on. The argument of the hypergeometric function in 
equation (jA.13p becomes 



A-y 



< 1 



7max(A+^) 

Thus, we can approximate 
A7 



This is the regime where the bulk of the electron 
population has already cooled and corresponds to the Fl4,p— l;p 
conditions 7 > 7b and 7 3> 7c or equivalently ^^-^ << 1, 



7n 



(A.18) 



(A.19) 
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The hypergeometric function with argument 1— — > This is the regime where the bulk of the electron population 
can also be approximated by unity. With the above approx- has already cooled and corresponds to the condition 7 > 7b 



imations, equation (jA.9[) reduces to: 
1 



or equivalently 7, > 7max- When 7 << 7n 

2 



cooled 



(A.20) 



T ^ r 
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Thus, the electron distribution in the cooled regime is given 
by the simple expression: 



p-1 y 2ao7 J 
Moreover the condition 7* > 7,nax leads to '^^ 



2ao7 



(A.30) 



< ^/r. 



cooled ■ 



ao(p~ 1) 



However, we can safely neglect the term 2^^, only if it is 
(A. 21) much smaller compared with the term y^. This condition , 
as discussed in section p.l.ip . is the same as 7 » 7c. Under 
this approximation: 



A.2. Magnetic field of the form B{r) = Bq^ 
The general solution of equation (|lip is 



rcoolcd 



7" 



-p+i 



iV(7,r) 



72 
7 



(A.22) N, 



cooled 



p-i 

aoip~ 1) 



r7 



-p-i 



(A.31) 
(A.32) 



or changing the integration variable 
1 
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A. 3. Constant magnetic field 
(A. 23) The solution of equation pT|) is: 

?,2 



r 7^ 

where constant ao is defined in section (|3.1.2|) . Assuming N{^,r) — / dfQ{j,r)-^ 
that the blast wave is in the decelerating phase (eq.([l])) the J~oo 7 

above equation becomes: 
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or changing the integration variable 
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where constant oq is defined in section (I3.1.3P . As described 
in the previous sections , the solution in terms of the inte- 
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If p ^ 2 the integral Ip can only be estimated in two 
regimes: 

— Uncooled 
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In this regime, 7 < 7c and 7* < 7max- Equation (jA.25|) 7* 
becomes 



1 2ao / 
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We outline next the points needed that will faciliate a 

co mparison of our analytic al results and those presented 

in lDermer fc Chiang! (|l998l) (h er eafter DC98). The expres- 

m , . . , . .1 . 1 o / ■ J. • 11 sions given by equations (IA.35I) . (IA.36I) coincide with these 

I akmg mto account that 7* > 7 and p > 2 (as it IS typically r /oo\ i-oa\ cT^no u 

■ ^- r \ A orr h 1 ,1 , of equations (33), (34) ot DC98. Ihis can be seen after 

taking into account the corresponding symbolism of our 

present work with the one used in DC98. The spatial co- 

(A.28) ordinates x,xq, used in DC98, and r, ro, used in our work, 

are identical. In the case of a spherical blast wave 



assumed), inspection of (jA.27p shows that 
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When the bulk of the electron distribution is considered to 
be relativistic then /? « 1 and 



P = /^7 ~ 7- 



(A.39) 
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The constant F used troughout our paper is related to 
To, xq of DC98 by 



(A.40) 



Our constant is related to the constant in the right hand 
side of Eqn. (33) of DC98 by the relation 



^c47rpoc 
nief 



(A.41) 



where = fc and po — mpUQ. The relationship between 
the two integrals Ipir) in DC98 and Ip{'y,r) in our paper 
is given by 



ao 



(A.42) 



Moreover, the constants vq and uj which appear in Eqn. 
(34) of DC98 can be expressed as 



vo = ao 

5fc 
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Also the variable r of DC98 can be expressed in terms of 
the distance r measured in the frame of the explosion by 
the relation 



T = — - r 



5cT 
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(A.45) 



Finally the exponent u which appears in the integral /p(t) 
takes the value — |. The characteristic Lorentz factor 7b 
which divides the analytic solution of equation (|A.35P into 
two branches is given by: 



7b 



2aQ [^5/2 _ y,5/2 

5cr ' 
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Taking into account 7max 3> 1 and r 3> equation (jA.46P 
becomes: 
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Therefore, in this case the two characteristic Lorentz factors 
of the distribution 7c, 7b coincide. In order to obtain the 
asymptotic solutions presented in section (|3.1.3p one follows 
the procedure presented in the previous sections. 
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